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Abstract 

We present an efficient generator for the process e~^e~ — > 4 fermions + 7 
through off-shell W pairs. It is based on a massless matrix element with 
leading 0{m'^) corrections. Only the resonant WW graphs are included. 
We have tested it against a matrix element without these approximations 
and found agreement to within ~1% at LEPII energies. 
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1 Introduction 



One of the main motivations of the LEPII electron-positron colhder is the study 
of pair-produced W bosons. In order for the anticipated high precision (~ 1%) of 
the measurements to be meaningful, they must be complemented by a complete 
calculation of the electroweak radiative corrections to one-loop level, which are of 
a similar order of magnitude. A subset of these are the hard radiative processes, 
in which a detectable photon is emitted. This process has been studied before for 
the production and decay of on-shell W bosons (for production see, e.g., [0, 0). 
However, the width of the W boson is such that one has to consider off-shell 
effects by looking at the complete process e~^e~ —>■ 4 fermions + 7. For most of 
phase space this is still dominated by the resonant W diagrams shown in Fig. |I[ 
The non-resonant W diagrams (an example is shown in Fig. ^) are suppressed 
by factors Tw/mw and can safely be neglected. The 77, 7Z and ZZ graphs 
shown in Fig. occur only when particle-anti-particle pairs are present in the 
final state. They also peak in a different region of phase space and can thus be 
excluded by cuts on resonant Z's and 7's decaying into observable particles. 

The non-collinear behaviour of the resonant graphs of Fig. |^ was already 
given in Ref. 0. We have extended this calculation to the full phase space 
(including mass effects |^) and all decay channels, added initial state coUinear 
bremsstrahlung and converted it into an efficient event generator with an 
interface to Jetset We also verified that the matrix element agrees with an 
independent approach p|, and checked that the contribution of the other graphs 
is indeed small. 

2 Method 

In this section we will sketch the way the program is built up and which methods 
are used to compute the various ingredients which make up the generator. This 
consists of the event generator package, a front end that chooses the final state 
particles, the phase space routines that transform the random variables to four- 
vectors and the matrix element, including mass effects and QCD corrections. 
Each of these parts is discussed in turn. Throughout we will use the notation of 
Fig. |. 

2.1 Event generator 

We use two approaches for the generation of weight one events which are useful 
for a full physics study.Q The first possibility is to throw random points into the 

^ There is also the possibility of using weighted events (normal Monte Carlo integration) for 
studies for which time-consuming tasks such as jet fragmentation and detector simulation are 
not needed. 
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Figure 1: The Feynman diagrams included in the generator. (The last row only 
occurs when the W decays hadronically) . 
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Figure 2: Some diagrams which are not included. Diagram a) is a non-resonant 
diagram, which is present in all channels but suppressed by a factor Tw/rriw 
Diagram b) occurs only in certain channels; it peaks in different regions of phase 
space than the W graphs. 




Figure 3: Assignment of the momenta 

integration region and accept an event if the ratio of its weight to the maximal 
weight exceeds yet another random number in (0, 1). This depends on a knowl- 
edge of the actual maximum weight, for which we take 20% more than the largest 
weight found in the first 500 events, which are discarded. As the integrand has 
been smoothed out by mappings so that the function does not fluctuate much 
this estimate normally exceeds the maximum weight. If a larger weight is found 
during the run the event is just accepted once for sure, and a second time with 
a probability depending on the ratio of the excess to the maximum weight. This 
approach has been implemented in the routine simplemc called by simple. It 
has an efficiency (events generated/events accepted) of about 5%. 

A more sophisticated method is to use the grid produced by an adaptive 
Monte Carlo algorithm such as Vegas or Bases P]. These adaptive algorithms 
use a discrete mapping on each integration variable to reduce the variance of 
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the integrand; generating points using the same discrete mapping improves the 
efficiency of the event generator. Event generation thus proceeds in two steps: 
first one finds a good grid by running the integration program, next one generates 
the events. The maximal weight has also been found already by the integration 
step. We integrated our programs with the Spring event generator, which builds 
on the Bases integration routine. The efficiency of this event generator can be as 
high as 20%, but this depends very much on the amount of time used to generate 
the grid. 

2.2 Front end 

The function which is called by the integration or event generation package is 
wwf . This takes as argument a point in a 13-dimensional unit hyper-cube and 
returns the weight associated with that point in picobarn. The corresponding 
four-vectors and particles are stored in common as described in section |4.3| for 
further treatment in the routine spevnt (jet fragmentation, detector simulation, 
etc.). 

The ffist activity is the conversion of the user-supplied values for the physics 
constants to coupling constants, W width and cut-off angles. This is repeated 
whenever these values change, printing the new values. Next the first semi- 
random variable is used to select the decay channel — leptonic, semi-leptonic or 
hadronic — from among the allowed channels specified by the user in wwf set 
(see section [4.1| ). After rescaling the variable, it is reused to select the leptonic 
decay channels (e, fi, r; again from those allowed by the user) and the hadronic 
channels (weighted according to the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements). One more choice needs to be made: the multi-channel integration 
algorithm means that the cross section is a sum of 4 (leptonic) to 6 (hadronic) 
different channels (in each channel the program focuses on the radiation off one 
of the external legs, see Ref. [§]). The same ffist random variable is again rescaled 
and used to choose the channel; finally the corresponding routine wgla-wg2b is 
called with 11 random variables. 

In the case that exponentiated collinear initial state bremsstrahlung is added 
to the full one-photon calculation the routine wwf ini computes the relevant struc- 
ture function factor and generates the initial state momenta after this radiation 
using the remaining random variable. All the momentum is assumed to be lost 
by just one of the incoming particles, although the function used includes up to 
two hard photons. This will only make a difference of order (with respect to 
the radiative cross section) when the cuts used discriminate between differently 
boosted events. 
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2.3 Phase space 

There are thus 6 different weight-computing routines, which each call their corre- 
sponding phase space routine f ivela-f ive2b. This phase space routine is called 
twice to construct a massive (pf) and massless (gf ) set of variables. The lat- 
ter is used to evaluate the matrix element consistently, the former to cancel the 
coUinear poles and to implement the experimental cuts. For the denominator, 
wgla-wg2b compute the jacobians of the other channels from the four- vectors 
with the routines gla-g2b. When the event satisfies some simple kinematical 
cuts in wwfrej the massless matrix element mat is computed. This is the same 
for all channels, but different mass corrections getiifflil-getiiim2 are added. Finally 
histograms are filled (wwf ill). 

The phase space routines evaluate the phase space element as follows. 

f ivela,gla: radiation off particle 1 (the incoming positron). The phase space 
integral is evaluated as 

dPS' = — ds^ — ds~ d PS2{s, , s~) d PS2{s^ jmljiril) d PS2{s~ ,ml,ml) , 

(2) 

where Ei, 9i and (pi are the energy and angles of particle i in the lab frame. 
In this equation s"^ are the momenta squared of the off-shell VF^ and 



dPS2{s,s+,s-)^d(P'^dcos0'^ ' ^^^^^^ ' , (3) 

with these angles defined in the CM frame. The Kallen function 

A is defined as A(x, y,z) = {x — y — zf' — ^yz. Analogous expressions hold 
for the other two 2-particle phase space elements. We map away the 1/t, 
l/E'r and — cos6'7) cx 1/ (piPy) behaviours of the matrix element, 

taking into account the cuts on the photon energy and angle in the first 
two. The mappings are defined in f ivela and mapt. 

fivelb.glb: radiation off particle 2 (the incoming electron). These routines are 
just the mirror image of the previous ones. 

f ive2a(0) ,g2a(0) : radiation off particle 4 (the outgoing e+, r+, J, s or h). 
The phase space integral is now evaluated as 

dPS2a = ^ds+ ^ds- dPS2{s,s+,s-) d PS^is^ , ml ml) dPS2{s-,mlml) , 

ZTT ZTT 

(4) 



5 



with the 3-particle Dahtz decay given by 



rl (,+ r»2 _ dE^ dE* d(t>lj ^^cicosg; 

with all quantities defined in the rest frame: E* is the energy of particle 
i, the angle between particles 3 and 7 (the photon), and (t>\,6*^ the 
orientation of particle 3 with respect to the W^. We now map away the 
1/t, l/E^ and l/(Vs+/2 -E* + (m^ -m2)/2Vs+) oc 1/(^4^7) behaviours of 
the matrix element, taking into account the cuts on the photon energy and 
angle as far as possible.0 Care was taken that the numerical stability is good 
enough for electron masses. The routines implementing these mappings are 
mapt, map el and mape2. 

f ive2a(l) ,g2a(l): radiation off particle 3 (the outgoing -u or c quark). These 
cases are the same as the previous ones with 3^4. 

f ive2b(0) ,g2b(0) : radiation off particle 5 (the outgoing e^, /i^, r^, s or b). 
These routines are obtained from f ive2a(0) ,g2a(0) by interchanging the 
W+ and W-. 

f ive2b(l) ,g2b(l) : radiation off particle 6 (the outgoing m or c anti-quark). 
These cases are the same as the previous ones with 5 6. 



2.4 Matrix element 

The massless matrix element was taken over from Refs 0, 0. It uses a massless 
helicity method, which breaks down in the coUinear region || pi 2. In this 
region we use a coUinear approximation described in the next subsection. The 
finite width of the W is introduced by replacing all propagators l/(p^ — m^r) 
by l/ijP' — + imy/^w) (without this the cross section is infinite). Q This is 
known to break gauge invariance, and is also not a unique prescription, as it 
depends on which factors have been cancelled against the numerator. The extent 
of the violation of electromagnetic gauge invariance had been checked before by 
changing the photon polarisation vector to its momentum. The gauge breaking 
terms turned out to be small: this test indeed gives a near-zero amplitude. 

We have checked this matrix element against one generated by the Grace 
package [Q] for all the possible final states.^ This matrix element includes masses, 
and gives us the option to test the contribution of non-resonant, ZZ, and 77 

^The directions of the momenta are unknown when the energies are chosen, so we have to 
err on the safe side; this gives rise to a loss of efficiency of about 4% at -^s — 190 GeV and 
slightly more at higher energies. 

•^The effect of using a p^-dependcnt width is negligible. 

*Due to the long evaluation time (more than one second on our workstation) this matrix 
element is not suitable for use in a Monte Carlo study, given the 13-dimensional phase space. 
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graphs. For massless particles and a zero W width the agreement was very good 
(5-8 decimal places). For a non-zero W width the two disagree by a small amount 
due to the non-unique way the width is introduced in each matrix element. In the 
non-collinear region the mass effects are also small, of order m^/ s as expected. In 
the Grace matrix element we also studied the effect of the gauge variant terms by 
explicitly varying the gauge parameters. The effect is significantly smaller than 
the accuracy of our calculation. 

The following tests were all performed at an energy y/s =190 GeV. For the 
integrated cross section with 0.1 GeV < < 60 GeV and no cuts on angles, the 
two methods agree to better than ~0.1%, for final states involving e, /i, u or d; 
better than ~0.25% for r's, and better than ~2% for b quarks. This deviation for 
b quarks occurs in the collinear region and will be hidden by jet fragmentation. 
With Grace we can also test the effect of the non-resonant graphs, which increases 
the above errors to <0.5% (leptonic without ZZ and Z7 graphs), <1% (udcsj) 
and <3% (udcb'j). The errors are somewhat larger in certain areas of phase space, 
but only when the amplitude is small. We have also tested that the exclusion of 
ZZ and Z7 diagrams in the case of, e.g., fiufiwy final states, is not a problem, 
given a small cut on such events coming observably from resonant Z's and 7's. 

2.5 Mass effects 

To introduce the correct collinear logarithms we multiply the massless matrix 
element by HAQiQi) / HAPiPy) , where the product runs over all charged particles. 
This, however, fails to take into account the double pole terms, which have a 
numerator proportional to but after integration give a finite contribution; 
these are added by hand. 

In the case of radiation off the initial state positron or electron (channels wgla 
or wglb) the matrix element is replaced in the collinear region {pik) < ml by the 
approximation (in getmml) 

|A.W|^-.4™„(il±f^-^)|A.<«)(MP. (6) 

with ^ = 1 — E-j / El the splitting variable and the matrix element without 
the extra photon. The mass term is added to the original matrix element in a 
much larger cone. The result agrees to within |% with the full massive matrix 
element 0. 

When the photon is predominantly radiated off the final state particle i (chan- 
nels wg2a and wg2b) we add to the matrix element the double pole terms (with 
getmm2) 

\M{s)\' - \M{s)\' -Ar^a,q^^^\M^\s)\' , (7) 

where qi is the charge of particle i in units of the electron charge. These terms 
are negligible outside the collinear region, and the sum agrees excellently (to 
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0{Tnf/s)) with the full massive matrix element, as can be seen from the tests 
cited in the previous section. 

2.6 QCD effects 

Leading order QCD effects have been incorporated in the rate by multiplying 
each hadronic W vertex by a factor (1 + Q;s(m^)/7r) with as(m^) = 0.133. 

3 Results 

We used the Vegas program to generate the following plots. The parameters 
used are mw = 80 GeV and T\y = 1.956 GeV; we work in the a scheme. No 
angular cuts have been applied and collinear initial state radiation is included. 
The photon energy spectrum in the leptonic channel at ^/s = 190 GeV is shown 
in Fig. |[ One would expect some structure around Tw as at this energy 

the W width starts to function as an infra-red cut-off; however, this is totally 
washed out by the radiation collinear to the light external charged particles. 

A plot of the total leptonic cross section with these parameters is given in 
Fig. |. To avoid the Z peak we took < E - (100 GeV)V(4£'). In the range 
of lower photon energy cuts shown, the cross section scales logarithmically with 
the cut-off, as it should. This only makes sense when this calculation is combined 
with the virtual and soft corrections, which will have the opposite logarithmic 
divergence to give rise to a finite total cross section. 

To give an idea of the attainable accuracy, we computed the cross section at 
each energy with 5 iterations of 10^ points; this gives a relative error of 0.2% and 
takes about 10 minutes per energy. 

4 Usage 

4.1 Parameters 

The physics parameters are divided into two classes. The constants which are 
unlikely to change are defined as parameters in the include file ww.h. These 
include the fermion masses, coupling constants and CKM matrix elements. The 
variable ones have to be defined in wwf set: these are the on-shell W mass (in 
GeV) mw and the physical width of the W wwidth. If wwidth is specified smaller 
than zero the one- loop on-shell width is computed. If schpar is zero the width is 
taken constant, otherwise an s-dependent expression is used. Finally one has the 
choice of working in the a-scheme (schpar=l) or Gj-scheme (2). This does not 
influence the final state photon couplings; this coupling constant is always taken 
to be a(0). 
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Figure 4: The photon energy spectrum at = 190 GeV in the leptonic channel. 
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Figure 5: The total leptonic cross section for > 0.1 GeV and E^>\ GeV. 
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The next group of parameters concern the decay channels which one wishes to 
be included. The parameter dectwo should be defined to be the sum of leptonic, 
semileptonic and/or hadronic. The leptonic decays of the and can 
further be specified with declpp and declpm respectively, where the possibilities 
are electron, muon and/or tau. There is no possibility to specify the hadronic 
decays. 

Finally one should specify whether one wants additional coUinear initial state 
photons included by setting inirad. 

4.2 Integration 

As the matrix element is infra-red divergent one has to specify (in wwf set) a lower 
limit on the photon energy of the generated events, ecut. Another problem occurs 
at the upper end of the photon energy range: if the energy of the photons (possibly 
including collinear initial state photons) is too large one reaches the Z peak. As 
the phase space routines do not take this peak into account the efficiency of the 
Monte Carlo will drop dramatically if this region is included. The contribution of 
the Z peak is negligible ((9(10~^), 1^), so this will not infiuence the final results. 
If eupcut is not defined in wwf set it is taken to be 60 GeV, which is reasonable 
for LEP II. A warning is printed if the Z peak is reachable. 

Often one is not interested in photons emitted collinear to final state particles 
or the beam pipe. To increase the generation efficiency one can implement these 
cuts in the phase space routines by setting acut (the angle between the photon 
and all charged particles) and bacut (the angle to the beam pipe) to a nonzero 
number of degrees. Note that the cone excluded around the beam pipe is the 
maximum of these two cuts. 

4.3 Events 

Whenever an event is accepted the routine spevnt is called. This stores the final 
state particles in jetset structures 0. This is also the place to include a further 
analysis or detector simulation. The four- vectors are stored in p (0:3, 7), the dot 
products in piDpj(7,7) (massive) and qiDqj(7,7) (massless) and the particle 
id's in iparti(7) in jetset codes. 

4.4 Summary of approximations 

As discussed before we made the following approximations in the event generator: 

- Only WW diagrams are included, no ZZ, 'jZ or 77 diagrams which could 
lead to some of the final states. With a small cut this is not a problem. 
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- Only the resonant WW diagrams are included; the effect of the non-resonant 
diagrams is minute in the case without extra photons P, |10] and we have 
tested it to be small here as well. 



- We use a massless matrix element, which has been corrected for collinear 
mass effects (the logarithmic and finite double pole terms). The neglected 
Oirrij/s) effects will be unobservable. 

- The W width is introduced in a gauge variant way; however, we have seen 
that the numerical effect of this is small in the LEP II region. 

- There should be an upper cut on the photon energy to avoid the Z peak; 
the numerical influence of this peak is completely negligible |p. 
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A Installation 



The generators have been written in reasonably standard fortran and have been 
tested on Sun (both SunOS 4 & 5), HP 700, DEC a and NeXT. The source code is 
available with anonymous ftp or WWW from pss058 .psi . ch in /pub/wwf . The 
included makefile should be easily adaptable to your site. It defines the targets 
wwfmc: integration with Vegas; wwfsimple: event generation with the simple 
algorithm; wwf bases: integration with Bases and wwf spring: event generation 
with Spring. Note that this distribution does not include a copy of Bases/Spring; 
these can be obtained from kekux.kek.jp in kek/minami/bases50. 
Some points to note are: 

- On OSF, the command ranlib does not exist; please delete these lines from 
the Makefiles. (HP-UX gives a laconic message.) 

- HP-UX does not have a flush routine used by the version of Vegas included; 
edit the lib/MaJsef ile to include the replacement flush. c. 

- As usual, the optimisation option on some of the compilers was found to 
be less than trustworthy. 
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